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Abstract 

We study the asymptotic behavior of a class of stochastic dynamics 
on interlacing particle configurations (also known as Gelfand-Tsetlin 
patterns). Examples of such dynamics include, in particular, a multi- 
layer extension of TASEP and particle dynamics related to the shuf- 
fling algorithm for domino tilings of the Aztec diamond. We prove 
that the process of reflected interlacing Brownian motions introduced 
by Warren in |Wj serves as a universal scaling limit for such dynamics. 

1 Introduction 

Consider N(N + l)/2 interlacing particles with integer coordinates xj, j = 
1, . . . , N, i = 1, . . . , j satisfying the inequalities 

x\_ x < x 3 r\ < x \. (l) 

and denote by GT^ the set of all vectors in Z N( - N+1 ^ 2 , which satisfy ([!]). 
An element of GT^ 5 -* is shown in Figure [l] Elements of GT^ are often called 
Gelfand-Tsetlin patterns of size N and under this name are widely used in 
representation-theoretic context. Whenever it does not lead to confusion, we 
use the notation x\ both for the location of a particle and the particle itself. 

In this article we study a class of Markov chains on GT^- 1 with simple 
block/push interaction between particles. Let us start with an example. Let 
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Figure 1: Interlacing particles. 



Y(t), t > be a continuous-time Markov chain taking values in GIT ^ and 
denned as follows. Each of the N(N + l)/2 particles has an independent ex- 
ponential clock of rate 1 (in other words, the times when clocks of particles 
ring are independent standard Poisson processes). If the clock of the particle 
yj r j n g S a t time t, we check whether Y?~ l {t— ) = Y?(t—)+l. If so, then noth- 
ing happens; otherwise we let Y/(t) = Yjj(t—) + 1. If Y/(t—) = Y^(t—) = 
■■■ = Y*g(t-), then we also set Y*g(t) = Y^(t-) + 1, . . . ,Y^(t) = 
Y£ k (t-) + 1. 

Informally, one can think of each particle having a weight depending on 
its vertical coordinate, with higher particles being lighter. The dynamics is 
defined in such a way that heavier particles push lighter ones and lighter ones 
are blocked by heavier ones in order to preserve the interlacement conditions. 

The Markov chain Y was introduced in |BF] as an example of a 2d growth 
model relating classical interacting particle systems and random surfaces aris- 
ing from dimers. The restriction of Y to the N leftmost particles x\, . . . , 
is the familiar totally asymmetric exclusion process (TASEP), while the par- 
ticle configuration Y(t) at a fixed time t can be identified with a lozenge tiling 
of a sector in the plane and with a stepped surface (see the introduction in 
[BF] for the details). 

More generally, we can start from any iV(iV + 1) /2 integer-valued random 
processes (in discrete or continuous time) with unit jumps of particles and 
construct a process similar to Y on through the same block/push 

interactions between particles (see Section [2] for a formal definition). We 
will refer to this process as the interlacing dynamics driven by the original 
iV(iV + l)/2 random processes. Some of such processes were studied in [BFJ, 
I WW] and (Nj. In particular, in the last article it is shown that interlacing 
dynamics driven by iV(iV-|- l)/2 independent simple random walks is closely 
related to the shuffling algorithm for the domino tilings of the Aztec diamond 
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studied in jEKLP] . 

In the present article we investigate the convergence to the scaling limits 
of the process Y and more general interlacing dynamics. The asymptotics of 
the fixed time distributions of such processes as t — > oo is well studied. Most 
results are easier to prove if the process is started from the densely packed 
initial conditions: Y/(0) = 1 — i, j = 1, . . . , N, i = 1, . . . , j. The reason is 
that such processes often have determinantal correlation functions (see |BFj ). 
However, many of the results hold in much greater generality. If one takes 
the limit t — > oo for a fixed JVeN, then the distribution of Y(t) converges 
to the so-called "GUE-minors" process, namely the joint distribution of the 
eigenvalues of an N x N Gaussian Hermitian matrix and the eigenvalues of 
its top-left corners (see [Bar] . |JN) . [OR]. [N]). 

If N goes to infinity together with t, then locally one unveils the Tracy- 
Widom distribution, the Airy point process and discrete point processes gov- 
erned by the sine kernel in the limit (see [BKj , [Jl] , [J2] ) , whereas the global 
fluctuations in this model are described by the Gaussian Free Field (see |BFj ). 
In the physics terminology the model fits into the anisotropic KPZ univer- 
sality class (see the introduction in |BFj and the references therein for more 
details). 

We perform the next step and study the joint distributions of Y and 
similar processes at various times. Currently we restrict ourselves to the case 
of fixed N. The limit process W- , j = 1, . . . , N, i = 1, . . . , j was introduced in 
jW] and can be constructed inductively as follows. is standard Brownian 
motion; given W- with j < k the coordinate is defined as a Brownian 
motion reflected by the trajectories W^Z\ (if i > 1) and W^~ l (if i < k) 
(see jW] or Section [2] for a more detailed definition). The process W has 
many interesting properties: for instance, its projection , i = 1,...,N 
is an iV-dimensional Dyson's Brownian motion, namely the process of iV 
independent Brownian motions conditioned not to collide with each other by 
means of Doob's /i-transform. 

Our main result is the following theorem. Let D([0, oo), ]R iV ( Af + 1 )/ 2 ) stand 
for the space of right-continuous functions from [0, oo) to M 7V ( Ar + 1 )/ 2 having 
left limits, endowed with the topology of uniform convergence on compact 
sets. 

Theorem 1. Let X(-; n), n = 1, 2, . . . be a sequence ofL N ( N+l ^ 2 -valued ran- 
dom processes with unit steps in each coordinate and let X{-] n), n = 1, 2, . . . 
be the interlacing dynamics driven by X(-;n), neN, respectively, with some 
initial conditions X(0;n) G GT^ ; ji6N, which are independent of the in- 
crements ofX(-] n), n £ N. Suppose that there exists a sequence {a n (-)} ng N of 
continuous real-valued functions on [0, oo) and a sequence {6 n }neN of positive 
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reals such that, as n — )■ oo, one has b n — > oo, the law of (X(-; n) — a n (-))/b n 
on D([0, oo),R N{N+1)/2 ) converges to the law of a standard N(N + l)/2 di- 
mensional Brownian motion, and the law of (X(0;n) — a n (0))/b n converges 
to some law W(0). 

Then the law of (X(-;n) - a n {-))/b n on D([0, oo), converges to 

the law of the process W with initial condition W(0) as n — >■ oo. 

Remark 2. If one wants to speak about discrete time processes, then it suffices 
to extrapolate them to continuous time processes with piecewise constant 
trajectories, to which one can apply Theorem [TJ 

Remark 3. Theorem [TJ for a special choice of driving processes X(-;n) was 
conjectured in jNj. 

Remark 4. Theorem [TJ is known to hold for certain projections of the process 
Y . Namely, in the literature one can find proofs of the convergence towards 
the joint distribution of Wf* , i = 1,...,N, the joint distribution of 
i = 1, . . . , N, the joint distribution of W{, i — 1, . . . , N and the distribution 
of W at any fixed moment of time. 

Using the well-known convergence of the standard Poisson process to the 
standard Brownian motion, we conclude from Theorem [TJ that the laws of 
the processes (Y(nt) — tn)/ 'y/n, t > converge in the limit n — > oo to the 
law of the process W, started from 0. Similarly, applying Theorem [T] to the 
interlacing dynamics of iV(iV + l)/2 independent simple random walks and 
using the convergence of the latter to standard Brownian motions, we prove 
the conjecture in [N] on the convergence of the particle dynamics related to 
the shuffling algorithm for the domino tilings of the Aztec diamond to the 
process W started from 0. 

The proof of Theorem [T] is based on the following idea. The limit process 
W can be constructed inductively using Skorokhod reflection maps in time- 
dependent intervals recently introduced in |BKR] , applied to the independent 
Brownian motions driving the particles. Moreover, the prelimit processes 
can be seen to be the images of the respective driving processes X(-;n) 
under similar Skorokhod maps. Putting these facts together with the joint 
continuity of the Skorokhod map in the driving path and the time-dependent 
boundaries, one obtains Theorem [TJ 

There are two directions in which Theorem [TJ might be generalized. First, 
it seems plausible that the condition on the prelimit processes to have unit 
steps in Theorem [TJ can be weakened. There are examples of the dynamics on 
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Gelfand-Tsetlin patterns in the literature (for example, interlacing dynam- 
ics driven by sums of geometric random variables as in [BFJ,|WWJ, or the 
shuffling algorithm for boxed plane partitions as in [BGJ, [BGRj, [Be]), for 
which the fixed time distribution is known to converge to the "GUE-minors 
process". Therefore, it is reasonable to expect that the laws of the paths of 
these discrete processes also converge to the law of W. One might also want 
to study what happens when the different components of the driving process 
have different asymptotic behavior, that is, when the sequences {a n (-)} n( z N , 
{b n }neN in Theorem [I] may also depend on the particles. Some particular 
results into this direction were obtained in Section 7 of |BG2] . 

In addition, Theorem [T] allows to recover many properties of the process 
W (for example, the ones studied in [W]) from the corresponding properties 
of the discrete prelimit processes, which are sometimes easier to prove (see 



e.g. |BF] ). For some more details in this direction, please see Remark 11 
below. 

It is worth noting that there is another natural way to construct a stochas- 
tic dynamics on which is based on the Robinson-Schensted-Knuth 
correspondence, see |01j . |Q2j . This dynamics is driven by N random walks 
as opposed to iV(iV + l)/2 in our case. While the evolution itself is very 
much different from the one considered in the present article due to strong 
correlations between components, the fixed time distributions (and also some 
other projections) are similar and can be made the same by an appropriate 
choice of the driving processes. 

The rest of the article is organized as follows. In Section [2] we introduce 
stochastic processes studied throughout the paper, in Section [3] we explain 
the relation of these processes to the Skorohod reflection map, and in Section 
[4] we prove Theorem [TJ 

Acknowledgements. The work on this article started while the authors 
were staying at MSRI and we would like to thank the organizers of the 
program "Random Spatial Processes". We are also grateful to A. Borodin 
and I. Corwin for useful discussions. V.G. was partially supported by RFBR- 
CNRS grants 10-01-93114 and 11-01-93105. 



2 Processes on interlacing particle configura- 
tions 

In this section, we give rigorous definitions of the stochastic dynamics on 
interlacing particle configurations that we study. To this end, let X(t) := 
{X\ (t), j = 1, . . . , N, i = 1, . . . , j}, t > be a random process taking values 
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m impose the following regularity conditions on X: 

1. The trajectories of X are cadlag, that is, for every t > 0, the limits 
X(t—) := lim s -|tX(s) and X(t+) := lim s ^X(s) exist, and X(t) = 
X(t+). 

2. Each coordinate X\ has unit increments, that is, for every t > 0, we 
have|X/(t)-X/(t-)|<l. 

Note that the regularity conditions imply, in particular, that on any time 
interval [0, T], the trajectory of X has finitely many points of discontinuity. 

Given a regular process X, we construct the interlacing dynamics X(t), 
t > driven by X and taking values in GT^* 1 as follows. The initial value 
X(0) can be any GT^ -valued random variable. If X is constant on some 
time interval [ti, t 2 ] , then so is X. If a t > is a point of discontinuity in 
X, then the value of X(t) depends only on X(t—) and X(t) — X(t—), and is 
given by the following sequential update. First, we define Xl(t), then Xf(t), 
i — 1,2, then X?(t), i — 1,2, 3, etc. To start with, we set 

x 1 1 (t) = x 1 1 (t-) + x±(t)-x 1 1 (t-y, 

in other words, the increments of the process X\ coincide with those of 
X\. Subsequently, once the values of X-(t), j = 1, . . . , k — 1, i = 1, . . . ,j 
are determined, we define Xf(t) for each i = 1, . . . , k independently by the 
following procedure, in which each step is performed only if the conditions 
of the previous steps are not satisfied. 

1. If i > 1 and X^(t—) = X*~^(t) - 1, then we say that particle X* is 
pushed by X^ and set X*(t) = Xf{t-) + 1. 

2. If % < k and X^t—) = X^ l (t), then we say that particle X t k is pushed 
by X*- 1 and set X t k (t) = X*(t-) - 1. 

3. If Xf (t) - Xf (*-) = 1, then we check whether i < k and ^ fc_1 (*) = 
X t k (t~) + 1. If so, then we say that particle Xf is blocked by A^ -1 and 
set X t k (t) = X^t-); otherwise, we set X t k (t) = X k {t-) + 1. 

4. If Xf (t) - X k (t-) = -1, then we check whether i > 1 and Xj^(t) = 
X k (t—). If so, then we say that particle X k is blocked by X k Z\ and set 
xf{t) = X k {t-); otherwise, we set X k {t) = X k {t-)-l. 

5. If X k {t) - X k {t-) = 0, then we set X k {t) = X t k (t-). 
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The above definition can be also adapted to discrete-time processes. If 
X(t), t = 0, 1, 2, . . . is such a process, then we define X(t), t > by 

X(t)=X(n) for n<t<n+l, n = 0,1,2,... (2) 

and let X(t), t = 0,1,2,... be the restriction of the interlacing dynamics 
driven by X to integer times. 

The following proposition enlists the properties of X which are immediate 
from the construction. 

Proposition 5. X(t), t > is a well-defined stochastic process taking values 
in GT^ . The paths of the process X depend only on the increments of the 
driving process X, but not on its initial value X(0). Suppose that X has 
independent increments, that is, for any T > 0, the process X(t) — X(T—), 
t > T is independent of the process X(t) — X(0), < t < T. Then X is a 
Markov process. 

Examples. 

1. If X is a family of N(N + l)/2 independent standard Poisson processes, 
then X is the Markov process Y defined in the introduction. 

2. If {£{(t) : j = 1, . . . , N, i = l,...,j, t = 0, 1, 2} is an array of i.i.d. 
Bernoulli random variables with parameter p and X{(t) = X/s=o£i(A)' 
then X gives the dynamics studied in [N] and [BF] . which is related 
to the shuffling algorithm for sampling random domino tilings of the 
Aztec diamond. 

Next we aim to introduce a continuous-space analogue of the above dis- 
crete processes. This process W was first defined and studied in [W]. It takes 
values in the continuum version of GT (A °, which we denote by GT { C N) . An 
element of is an array of reals x{, j = 1, . . . , N, i = 1, . . . ,j satisfying 

the inequalities 

< xlZl < x\. (3) 

The process W taking values in GT^ is defined level by level. For 
j = 1, we let Wl be a standard Brownian motion. For j > 1, we define the 
component W- as a standard Brownian motion (independent of those used 
in the previous steps) reflected on the trajectories of W?Zi and W?~ . In 
particular, Wj and W( are reflected on only one trajectory. We can start W 
from any deterministic or random initial condition with the only restriction 
being that W(0) belongs to GT^ with probability 1. For more details, 
please see [W], Section 2] and (WJ Section 4]. 
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Properties. As shown in [W], the process W started from possesses certain 
interesting properties. In particular, the projection W t N , i = 1, . . . , N evolves 
according to the Dyson's Brownian motion of dimension N. Also, for any 
fixed t, the distribution of W it) is given by the distribution of the eigenvalues 
of corners of an N x N Gaussian Hermitian matrix (see |JNj . |OR] ) . 

Another intriguing property is the following Markovity: Fix an iV G N 
and < t\ < t 2 < ■ ■ ■ < tjy and consider the process Z which is defined 
to coincide with the Revalued process (W* : % = 1, . . . , k) on the time 
interval [tk-i,tk) (here, we use the convention to = 0). Then the process 
Z{t), < t < is Markovian and its transitional probabilities can be given 
explicitly, see [W] proof of Proposition 6] and also |BFl Proposition 2.5] for 
a similar statement in discrete settings. 

3 Skorokhod maps 

The proof of our main result relies on the observation that the limiting pro- 
cess W is given by the image of the vector of driving Brownian motions, 
started in W(0), under an appropriate Skorokhod reflection map, whereas 
the prelimit processes X(-;n), n G N are given by images of the driving pro- 
cesses X(-;n), n G N, started in X(0;n), n G N, respectively, under slightly 
modified Skorokhod maps. Our construction relies on [BKRj Theorem 2.6], 
which we state in the following proposition for the convenience of the reader. 

Proposition 6. Let I, r be two right- continuous functions with left limits on 
[0, oo) taking values in [—00,00) and (—00,00], respectively. Suppose further 
that l(t) < r(t), t > 0. Then, for every right- continuous function ip with left 
limits on [0, 00) taking values in M, there exists a unique pair of functions 4>, 
rj in the same space satisfying 

1. For every t G [0, 00), (p(t) = ip(t) + rj(t) G [l(t), r(t)]; 

2. For every < s < t < 00, 



T)(t) — T](s) > if ' 4>{u) < r(u) for all u G (s, t]. 
rf(t) - i](s) < if(p{u) > l(u) for allu G (s,t\; 



(4) 
(5) 



3. For every < t < 00 



r){t)- V (t-)>0 if<j>(t)<r{t), 
r)(t)-r)(t-)<0 if (Pit) > lit), 



(6) 
(7) 



where r)(0— ) is to be interpreted as 0. 
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Moreover, the map T, which sends the triple (l,r,i/j) to the function <fi is 
jointly continuous with respect to the topology of uniform convergence on 
compact sets. 

We will refer to T(/,r, ■) as the extended Skorokhod map in the time- 
dependent interval [l(t),r(t)], t > 0. It will not be important in the following, 
but it is worth mentioning that the map T can be given explicitly (see |BKP4 
Theorem 2.6]). 

3.1 Skorokhod maps for the prelimit processes 

Since the definition of X(-;n) does not depend on the initial value X(0;n) 
we may assume that X(0;n) = X(0;n). In this subsection we show that in 
this case each of the prelimit processes X(-;n), n G N can be constructed 
as the image of the vector of driving processes under a deterministic map, 
which we will refer to as the discrete Skorokhod map. 

In our construction we fix an n G N and proceed by induction over the 
number of levels N. For N — 1, we set X(-;n) = X(-;n) = X(-;n). For 
N > 2, we may assume that the process X(-;n) with (N — 1) levels, with 
initial condition being the restriction of ^(0; n) to the first (N— 1) levels, has 
already been constructed. We fix a path of this process and will construct 
the corresponding paths of the particles on level N . 

We consider first an 1 < i < N and define X t N (-;n) as the image of 
X^(-;n), started in Af/^O; n), under the extended Skorokhod map in the 
time-dependent interval 



■;n) as the image of X^(-;n), started in X^(0;n), under the extended 
Skorokhod map in the time-dependent interval 



in the sense of Proposition |6j Similarly, for i = N, we define n) as the 

image of X^(-;n), started in X$(0;n), under the extended Skorokhod map 
in the time-dependent interval 




(-oo.A^-^n) -1], t>0 



(9) 



[X^(t;n),oo), t>0 
in the sense of Proposition |6j This finishes the construction. 



(10) 
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We remark at this point that the paths of the processes X(-;n), n G N are 
given by images of the paths of the driving processes X(-;n), n G N, started 
in X(0;n), n G N, respectively, under a deterministic map $ depending 
only on N (and not on n): 



X(.-n) = $ 6K (X(-;n)), nGN. (11) 

We will refer to Q SK as the discrete Skorokhod map. Finally, for each nGN, 
we define y(-',n) as the process constructed from x ('' n )~ a "(-) by employing 
the discrete Skorokhod map in the rescaled coordinates; that is, repeating 
the procedure above, but rescaling all initial conditions, processes involved 
and the constant 1 in ([8]), (|9| according to the change of coordinates 

[0, oo) x Z N ^' 2 -> [0, oo) x R^+W 2 , ( t> x ) h. ^ x ~^(t) \ (12) 

The following proposition shows that Af(-; n), n G N are in fact the processes 
defined in Theorem [lj and that y(-', n), n G N coincide with the correspond- 
ing rescaled processes. 

Proposition 7. Suppose that the processes X(-\ n), n G N and X(-; n), n G N 

are constructed by using the same driving processes X(-;n), nGN. Then, 



X(-;n) = X(-;n), nGN, (13) 



y(-;n) = X ^ n) r^'\ neN (14) 



with probability 1. 



Proof. Fix an n G N. We start with the proof of (13) for that n and proceed 



by induction over the number of levels N. For N — 1, the statement (13) 



follows directly from the definition of X(-;n). Now, let N > 2 and suppose 



that (13) holds for all N' < N. Then, the paths of the particles on the first 
(N — 1) levels of X(-] n) must coincide with the paths of the particles on the 
first (N — 1) levels of X(-; n). Now, consider the trajectories of the particles 
on the iV-th level in the two processes. If the trajectories coincide up to some 
time t%, then they clearly coincide up to time t 2 > t\ as long as no particles 
are pushed or blocked in X(-; n) in the time interval [ti,^]- 

In case that a particle on level N in X(-;n) is pushed, the trajectory of 
the pushing particle must have a jump of size 1 at the moment of the push 
and the other particle is pushed accordingly. The definitions of § SK and 
of the extended Skorokhod map in a time-dependent interval incorporated 
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in Q SK show that the value of X(-;n) = Q SK (X(-;n)) after the push must 
coincide with the value of X(-;n) after the push. 

Similarly, when a particle on level N in X(-;n) is blocked, the blocking 
particle must be at distance 1 from the blocked particle at the moment, 
when it blocks the other particle, and the trajectory of the component in 
X(-; n) driving the blocked particle must have a jump of size 1 at the same 
moment. By the definitions of $ and the extended Skorokhod map in a 
time-dependent interval, the value of X(-;n) = § SK (X(-;n)) after a particle 
is blocked must coincide with the value of X(-;n) at the same moment of 
time. 



The statement ( 14 ) can be shown by arguing as in the proof of the state- 



ment ( 13 ), but working in the coordinates obtained by the rescaling ( 12 ). □ 



3.2 Skorokhod map for the limiting process 

In this subsection we will employ Proposition [6] to show that the process W 
can be viewed as the image of the vector of driving Brownian motions, started 
in Jy(0), under a deterministic map, which we will refer to as the continuum 
Skorokhod map. To this end, we let B := {B^, j = 1, . . . , N, i = 1, . . . , j} be 
the collection of independent standard Brownian motions driving the process 
W and will construct an auxilliary GT^-valued process W on the same 
probability space, which a posteriori will turn out to coincide with W . 

The construction proceeds by induction over the number of levels N . For 
N = 1, we set W = B\ = W . For N > 2, we may assume that the paths of 
the particles on the first (N — l) levels of the process W, with initial condition 
being the restriction of W(0) to the first (N — 1) levels, have already been 
constructed. Then, for 1 < i < N, we define the path of the particle as 
the image of the driving Brownian motion Bf, started in Wf^lO), under the 
extended Skorokhod map in the time-dependent interval 

[W^\t\Wr\t)l t>0 (15) 

in the sense of Proposition [6j For i — 1, we define the path of the particle 
Wi as the image of the driving Brownian motion Bf , started in ^^(0), 
under the extended Skorokhod map in the time-dependent interval 

{- oo,W N -\tj\, t>0 (16) 

defined as in Proposition [6] Finally, for i = N, we define the path of the 
particle W$ as the image of the driving Brownian motion B^, started in 
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Wjy (0), under the extended Skorokhod map in the time-dependent interval 

[W^(t], oo), t>0 (17) 

in the sense of Proposition |6j 

We note that the process W is given by the image of the Brownian motion 
B, started in W^O), under a deterministic map ty SK : 

W = ^ SK (B). (18) 

We will refer to ty SK as the continuum Skorokhod map. The following propo- 
sition shows that W coincides with W. 

Proposition 8. Let the processes W and W be defined on the same proba- 
bility space and be driven by the same collection B of independent standard 
Brownian motions. Then, with probability 1: 

W{t) = W[t) for all t > 0. (19) 

Proof. We proceed by induction over the number of levels N. For N = 1, the 
statement of the proposition follows directly from the definition of W. Now, 
consider an N > 2 and assume that the proposition holds for all N' < N. 
Since the paths of the particles on the first (N—l) levels in W and W coincide 
with the paths of the particles in the respective processes with (N — l) levels, 
we have W- = W- for all j = 1, . . . , iV — 1, i = 1, . . . , j with probability 1. 

Now, fix an 1 < i < N and consider the particle W t N on the iV-th level in 
W and the corresponding particle W t N in W. Since both W t N and W t N solve 
the extended Skorokhod problem in the time-dependent interval 

[W^^W?- 1 ®] = [WiLffaW?- 1 ®], t > (20) 

for the Brownian motion B-, started in 14^(0), by their respective construc- 
tions, the uniqueness statement of Proposition [6] implies that it must hold 
W t N = W t N with probability 1. In the cases that % = 1 or i = N, one can 
argue in the same manner, invoking the uniqueness statement of Proposition 
[6] for the extended Skorokhod maps in the time-dependent intervals 

(-oo,W 1 N - 1 (t)] = (-oo,W 1 N - 1 (t)}, t>Q, (21) 
[W»:t(t),oc) = [W%lt(t),oo), t > 0, (22) 

respectively, to conclude that Wf = W± and W$ = W$ almost surely. 
This finishes the proof of the proposition. □ 
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Remark 9. Although this will not play a role in the following, we remark that 
due to Proposition [8] and the findings in [W] the particles on the same level in 
W never collide after time 0. Thus, if at time no two of these particles are at 
the same location, |BKP4 Corollary 2.4] shows that the extended Skorokhod 
maps used in the construction of W are in fact regular Skorokhod maps in 
the sense of p<TTl Definition 2.1]. 



4 Limiting procedure 

We now come to the proof of the main result of the paper, stating that the 
process W serves as a universal limiting object for discrete processes with 
block/push interactions. 

Proof of Theorem [7| We proceed by induction over the number of levels 
N. For N = 1, there is nothing to show. We now let iV > 2 be fixed and 
assume that the theorem holds for all N' < N. Since the first (N — 1) levels 
of the processes M^lzMi ; n £ evolve as the corresponding processes 
with (N — 1) levels, we conclude that their laws converge to the law of 
{W? : j = 1, . . . , iV — 1, i = 1, ... ,j} on D([0, oo), M^W 2 ), started 
according to the restriction of W(0) to the first (N — 1) levels. By the 
Skorokhod Embedding Theorem in the form of |Dul Theorem 3.5.1], one 
can define the processes on the first (N — 1) levels of X{- \ n) together with 
the processes Xf(-;n), i = 1,...,N, started in X t N (0;n), i = 1,...,N, 
respectively, for all values of n 6 N on the same probability space such that 

a n (-) n^ wl J = i,...,iV_i, , = 1,... ;J , (23) 
Xf(-;iV)-a„(-) n^oo^jv i = 1> _ >JV (24 ) 

On 

almost surely, where \W{ : j = 1, . . . , N — 1, % = 1, ... ,j} is a GT^ -1 ^- 
valued process of the law described in Theorem [l] and (B^, . . . , B§) is an 
N- dimensional standard Brownian motion started in (W^(0), . . . , W$(0)). 



Now, recall from Proposition [7] and the definition of the processes y(-;n) 



n G N that the paths of the processes — ^'"^ } i = 1, . . . : N, n & N are 



given by the images of the driving processes Xi - - , i — 1, . . . , N, n E 

N under the extended Skorokhod maps in the appropriate time-dependent 
intervals (we will refer to the latter as If{n)). Moreover, from Proposition [8] 
and the definition of the process W we know that the paths of the processes 
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W { , i = 1,...,N are given by the images of the driving Brownian motions 
Bf, . . . , B^j under the extended Skorokhod maps in the appropriate time- 
dependent intervals (we will refer to the latter as If (pa)). 

It now suffices to observe that, for all 1 < i < N, the left (resp. right) 
boundaries of the time-dependent intervals lf{n) converge in the limit as 
n — > oo to the left (resp. right) boundaries of the time-dependent intervals 
lf(oo) in D([0, oo), M) with probability 1. We remark that here the assump- 
tion b n — > oo as n — > oo is used. By applying the continuity statement 
in Proposition j^] we conclude that the paths of the processes — - an<y - , 
i = 1, . . . , N converge in the limit as n — > oo to the paths of the processes 
Wf , i = 1, . . . , N almost surely in D([0, oo), R). This finishes the proof. □ 

Remark 10. The proof of Theorem [T] shows that the result of the theo- 
rem is universal in the following sense. Suppose that X(-;n), n G N are 
such that the processes x ('' n )~ a "(') n E N converge in law to an arbitrary 
D([0, oo), ]R 7V ^ Ar+1 ^ 2 )-valued process B as n — > oo, where {a n (-)} n£ N is a 
sequence of continuous real-valued functions on [0, oo) and {&„}„ e N is a se- 
quence of positive reals such that b n — » oo as n — > oo. Then, the laws of the 
processes x ('> n )~ an t) n (zN converge to the law of the image of B under the 

On 

continuum Skorokhod map ^ SK in -D([0, oo), M JV ( Ar + 1 )/ 2 ) as n — > oo. 

Remark 11. Theorem [T] shows that one can reprove the properties of the 
process W and its projections stated at the end of Section [2] by using the 
corresponding properties of the process Y and taking limits. Indeed, one 
can reexpress the transition probabilities of the process W as limits of the 
corresponding transition probabilities of the appropriately rescaled versions 
of the process Y. 
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